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The multiscale entanglement renormalization ansatz is applied to the study of boundary critical 
phenomena. We compute averages of local operators as a function of the distance from the boundary 
and the surface contribution to the ground state energy. Furthermore, assuming a uniform tensor 
structure, we show that the multiscale entanglement renormalization ansatz implies an exact relation 
between bulk and boundary critical exponents known to exist for boundary critical systems. 



PACS numbers: 03.67.-a,05.30.-d,89.70.-a 

I. INTRODUCTION 

Variational approaches based on tensor networks pQ 
are a novel powerful numerical tool believed to be the key 
ingredient to simulate efficiently quantum-many body 
systems. Although a detailed understanding of their 
potentialities and their limitations is presently under 
scrutiny, there are already a number of encouraging re- 
sults. Variational Tensor Networks (VTN) are free of 
most of the problems of traditional numerical methods. 
Differently from quantum Monte Carlo methods, VTN 
do not suffer of the sign problem. Compared to density 
matrix renormalization group [2], they are more versa- 
tile and allow to simulate efficiently critical correlations, 
long-range interactions and two- and higher-dimensional 
quantum systems. Indeed the density matrix renormal- 
ization group can be reformulated in terms of a partic- 
ular class of tensor networks known as Matrix Product 
States [3 . VTN include also projected entangled pair 
states [4] that generalize matrix product states in dimen- 
sions higher than one and weighted graph states [5] de- 
signed to study systems with long-range interactions. 

Among the proposed VTNs, a very promising one is 
the so-called Multiscale Entanglement Renormalization 
Ansatz (MERA) [6]. MERA has been already applied 
successfully to the study of a number of different phys- 
ically relevant systems, like quantum models on a two- 
dimensional lattices [TJ |8] , interacting fermions [9 , and 
critical systems [TQl [TTJ [I2j [13], only to mention a few 
of the most remarkable examples. The capability of 
MERA to describe accurately critical systems derives di- 
rectly from the scale-invariant self similarity of its tensor 
structure, intimately related to a real- space renormal- 
ization procedure. The structure of the MERA state is 
designed [6 in such a way to reproduce scale-invariance 
and so, in one-dimensional systems, it naturally encoded 
several important features of the Conformal Field The- 
ory (CFT) underlying the critical lattice model [11 . The 
critical exponents can be computed directly from the 
spectrum of the MERA transfer matrix [TP] . 



Critical systems can however lack translational invari- 
ance, due to the presence of a physical boundary or to 
an impurity. The study of boundary critical phenomena 
is, since many years, a very active field of research which 
ranges from the study of the critical magnets with sur- 
faces to quantum impurity problems (as e.g. Kondo) or 
the Casimir effects (for a review of the field see for exam- 
ple [H]). The presence of the boundary does not spoil 
conformal invariance. Oppositely, boundary CFTs have 
a very rich structure and a deep mathematical foundation 
(see e.g. [33]). While in fact it started only as the study of 
critical two-dimensional systems in system with bound- 
aries (surface critical behavior), it found applications to 
open-string theory (D-branes), quantum impurity prob- 
lems [16 , quantum out-of-equilibrium studies (quantum 
quenches [T7] ) just to cite a few. Furthermore, it at- 
tracted a large attention from the mathematical commu- 
nity for the recent developments of stochastic Loewner 
evolution [T5] . 

In view of the connection between MERA and CFT, it 
is natural to wonder whether the MERA tensor-network 
can be employed to study quantum systems with bound- 
aries. In this work we introduce and analyze an entangle- 
ment renormalization tensor-network design which takes 
into account the presence of a (critical) boundary, and 
we study its properties. This is implemented by allowing 
the edge spin at the boundary of the system to interact 
at each step with an ancillary element, describing a fic- 
titious degree of freedom. Similarly to Refs. [10J [TTJ [12] , 
we will focus on homogeneous configurations, where ten- 
sor elements of the same class are also identical to each 
other. Interestingly enough this ansatz is able to capture 
some important properties predicted by boundary CFT. 
Specifically, we show that the critical exponent associated 
to the decay of any one-point function (as function of the 
distance from the boundary) is always half of the one of 
the bulk two-point correlation function corresponding to 
the same scaling operator. We also compute the bound- 
ary corrections to the ground state energy. 

The paper is organized as follows. In Sec. [II] we in- 
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Figure 1: Left panel - Entanglement renormalization network representation for N = 32 sites. In the bulk the MERA structure 
consists of isometries and disent anglers. The boundary is represented by an additional ancillary system indicated by the grey 
stripe. Right panel - Alternative insertion of the boundary obtained by applying the ancilla interaction tensors (magenta 
crescent-moons) at the same level of the disent anglers instead of the renormalizers, as done in the left panel. 



troduce the tensor network and its main properties. In 



Sec. [Ill] we discuss how the expectation values of local 
observable can be computed. Assuming a uniform tensor 
structure, in the thermodynamic limit, these expectation 
values decay as power law. We relate the associated crit- 
ical exponents to those of the corresponding two-point 
correlations in the bulk. In Sec.[lV|we discuss the bound- 
ary corrections to the ground state energy. The conclu- 
sions of our work are summarized in Sec. [V] 



II. THE TENSOR NETWORK 

Consider a ID lattice of N = 2 n+2 spins (sites), of a 
given local dimension <i, with open boundary conditions. 
A generic pure state of such system can be expressed as 



(1) 



state (yellow squares) whose bonding dimensions coin- 
cide with the coordinate space of the ancillas, which is 
coupled to the bulk via local coupling-elements (drawn as 
magenta crescent-moons in the figure). Via purification, 
we can always choose the dimension of such ancilla to be 
large enough so that the resulting interaction is described 
by a unitary operator, that we indicate as a, a type-( \ ) 
tensor of elements a^ 1 '^ 2 . Similarly to the case of the 
As and of the x s > we will also assume these elements to 
be uniform in the network (possibly allowing the ones on 
the left-hand-side of the structure to differ from the ones 
on the right-hand-side [19]). 

As customary with MERA-like configurations, to en- 
force efficient evaluation of local observables and correla- 
tion functions, the various elements of the network are 
assumed to obey specific contraction rules (a detailed 
analysis of the efficiency requirements for MERA can be 
found in Refs. [6j [I2j [20]). In particular the renormal- 
izers and the disentanglers obey isometric and unitary 
constraints respectively, i.e. 



with a canonical basis for the single qudit and T 

a type-( ° ) tensor. Following the prescriptions of the 
MERA structure [6 , we assume a formal decomposition 
of T which is schematically sketched in Fig. [I] Here we 
use the standard graphical convention (see e.g. Ref. [6]) 
for which each node of the graph represents a tensor (the 
emerging legs of the node being its indices), while a link 
connecting any two nodes represents the contraction of 
the corresponding indices. The yellow element on the top 
of the figure describes a type-( ° ) tensor C of elements 
C a ,^i^2,"- that we can call hat tensor. The green 
triangles represent instead the same d x d 2 renormalizer 
tensor A of type-( ^) of elements A^ 2 , and the blue 
circular elements represent the same d 2 x d 2 disentangler 
tensor \ of type-( \ ) of elements Xe^e^ • 

At the boundary, we introduce extra tensors that cou- 
ple the sites at the border with an ancillary degree of free- 
dom represented in Fig. [l]by the thick grey strip. These 
new elements form the lateral edges of the network and 
describe the boundary at each level of the MERA, i.e. at 
each level of the renormalization flow. As shown in the 
figure, each of them can be viewed as a matrix product 



E\u \ki,k 2 
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k\ ,k 2 



ki ,k 2 



where Sf is the Kronecker delta, while A^ 



(\£ \* 

\ U\ ,u 2 J 



and Xe^e^ = (Xui'fn 2 )* are tne adjoint counterparts of 
the A and x tensors respectively, obtained by exchanging 
their lower and upper indices and taking the complex 
conjugate. Similar conditions are imposed also for the 
edge tensors 



k\ ,k 2 



ui ,u 2 -ki 
a k 1: k 2 a h 



= S?} si 



(3) 



These rules are graphically represented in Fig. [2] 
Finally to ensure proper state normalization, 
the tensor C is supposed to satisfy the identity 

EC r i ,c a ^^'" 5^4, a' _ 1 

It is worth noticing that, by simply re- arranging the 
various tensorial components, an entanglement renormal- 
ization configuration which differs from the one given in 
Fig. [l] can be obtained. In fact, the ancilla interaction 
tensors a (the magenta crescents) can be applied at the 
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Figure 2: Rules for the elements of the tensor network of 
Fig. [I] The black marks drawn upon the tensors are used 
to make distinction between top and bottom. The first two 
schemes correspond to the the contractions of Eq. Q, while 
the last one represents the unitary constraint for the edge 
elements of the tensor. 

same level of the disent anglers, instead of the renormal- 
izers (see Fig. [I] right panel). Despite their different ap- 
pearance, it can be shown that these two structures are 
formally equivalent to each other. This can be verified 
by grouping together the edge-ancilla interaction with the 
nearest linked element belonging to the lower half-level. 
Now, by just performing a polar decomposition [2Tj . we 
obtain a structure of the opposite type (although, the 
very first spin of the chain is now taken out the system 
and put into the ancilla, while the second one becomes 
the edge spin). Having acknowledged this equivalence, in 
the rest of the paper we will work with the structure of 
the left panel of Fig. [I] 

III. LOCAL AVERAGES IN THE PRESENCE 
OF A BOUNDARY 

In the presence of the boundary, the average of any lo- 
cal observable depends on the distance from the bound- 
ary itself. For a one-dimensional critical system, the case 
we consider here, the space-dependence will be a power- 
law characterized by a set of critical exponents. In this 
Section we show how to compute local averages and how 
to extract these exponents. 

Consider a family T = {|?/^ n )); n = 2, 3, • • • } of states 
|^ (n) ) of increasing sizes, described by homogeneous net- 
works of the form shown in Fig. [TJ each sharing the same 
structural elements (renormalizer, disentangler, edge- 
ancilla interaction, hat). For such family we want to 
calculate the expectation value of a general observable 
acting on a small group of neighboring sites located at a 
given distance from the closest edge of the system, say 
the left one. For instance, in the case of a three-sited 
observable 6^ acting on the sites £,£ + 1, and i + 2 [22] 
we have 

<e<><»> = <v.<">|e*iv (n) > = Tr[e t ■ p# +m+2 ], (4) 

where the site indices are counted starting from the left- 
most spin as the first one, and where pfj+i is the 
reduced density matrix of \^ n ^) associated with the se- 
lected spins. 

We assume a uniform MERA structure. This assump- 
tion may seem an over-simplification for a system which 



is not translational invariant, but it turns out that it 
naturally accounts for the underlying (boundary) con- 
formal invariance. From the locality requirements im- 
posed in Fig. |2j it is straightforward to verify that for all 
1 ^ i ^ 2 n+1 — 2 and n ^ 1, the following recursion rules 
apply: 

O) —Ti( ( n-1 ) ^ 

P 21, 21+1,21+2 — U L \ Pl,l+\£+2 J ' 

(5) 

M — T) { ( n_1 ) ^ 

^2^+1,2^+2,2^+3 — Ur \ P 1,1+1,1+2 J ' 

where p^+1^+2 anc ^ P^e+i e+2 are 3-sites reduced density 
matrices of \^ n ~^) E T . In these expressions, Vl and 
Vr are completely positive trace preserving (CPT) maps 
that depend only on the bulk elements of the network (in- 
deed they coincide with the Vl and Vr maps of an ordi- 
nary homogenous MERA with the same A and x [ID]) and 
whose formal expression is graphically depicted in Fig. [3] 
By means of the the renormalization procedure implied 
by Eq. ([5]), at each application of the map the site over 
which the average is performed approaches the boundary 
in an exponential fashion. Correspondingly the network 
depth decreases linearly. Upon reaching the boundary 
one has to define further operations: 

Pl?2,3 = (P%1$) > ( 6 ) 

where A refers to the degree of freedom belonging to the 
ancillary system, and JCl is again a CPT map, sketched 
in Fig. [4] (left). At this point the causal cone of the ancilla 
jointed with the first two sites is stable. Indeed one has 

& = Bl (ptJ) - (7) 

where Bl is the CPT map represented in Fig. [4] (right). 
Analogously to Bl and JCl, we define the CPT maps at 
the right boundary from the mirror images of Fig. [4] and 
we call them Br and JCr respectively. 

Because of the stability of the causal cone [6j [TSJ [20] , 
approaching the thermodynamical limit we can deter- 
mine the reduced density matrix in proximity of the 
boundary by calculating the fixed point of Bl- This is 
unique provided that the CPT map is mixing [TQl [23] . 
i.e. 

lim p { 2\ , 2 = p f A ^ 2 = B L (pa,i, 2 ) • (8) 




Figure 3: Sketch of the CPT maps a) Vl and b) Vr; the 
contracted tensors picture represents their Kraus decomposi- 
tion. 
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Figure 4: Tensor graph representations of the CPT map JCl 
(left) and Bl (right) defined in Eq. (J6| and |7| respectively. 



We can now use this argument to obtain the expectation 
value in Eq. Q for infinite volume. The resulting expres- 
sion becomes particularly simple when Vl = Vr = V. In 
this case indicating the integer part of log 2 £ with |_log 2 £\ , 
we have 



Tirfe-I>Liog 2 <J o)C L (p f A12 ) 



(9) 



where 



stands for super-operator composition and 



where V T describe a r reiterated applications of the map 
V. 

We can now exploit the Jordan block decomposi- 
tion [21 to simplify further this expression. Adapting 
the derivation for the bulk in Ref. [12] to the boundary 
case, we easily get 



(CO) 



(10) 



where the sum spans over the eigenvalues noiV and g K (-) 
is a (finite degree) polynomial in its main argument with 
coefficients which depends on p^ A12 an d T>. Since CPT 



maps are contractive, the values of n entering in Eq. (10) 
belong to the unit circle (i.e. \n\ ^ 1). Furthermore, 
if V is mixing (which is a reasonable assumption [TO] ) 
then its spectrum admits a unique eigenvector (the fix 
point pf of the channel) associated with k = 1; all the 
remaining eigenvalues have modulus strictly smaller than 
1. Under these circumstances, in the limit of large dis- 
tance from the boundary, the quantity (B^ 00 ) converges 
toward its bulk limit which is obtained by computing the 
expectation value of O on the fix point of the channel, 
i-e. (Obuik)^ = Tr[6p/]. The deviations from such lim- 
iting expression can be evaluated by keeping the largest 
contribution associated with the terms with n ^ 1. This 
yields an exponential decay in log 2 £ of the form 



(e bulk )(°°) + ^ K Liog^j ffK ( L log 2 £J) 

(©bulk 



,(oo) 



»"V(log 2 *), (11) 



where k is the eigenvalue of V which has the largest ab- 
solute value smaller than one and which contributes not 



trivially to Eq. (10), and where g'(\og 2 £) is instead some 



complicated function which is dominated by a polyno- 
mial of log 2 £. In particular, if O is an eigenvector of the 
Heisenberg adjoint of the channel P, then Eq. (11) yields 



an exact power-law decay (without polynomial correc- 
tions), i.e. 



(O 



p>„\(°°) 



Jlog 2 «J 



(12) 



where c = Tr[6 K>l(Pa i 2)]' an ^ wnere ^ is the as- 
sociated eigenvalue (notice that for such a O one has 
(©buik)^ = 0). The above expressions show that the 
quantities — log 2 \k\ play the role of the critical exponents 
of the system. It is now evident that these critical expo- 
nents are the half of the corresponding ones for two-point 
correlation functions computed in the bulk, a well-known 
result in conformal field theory [15]. For instance fixing 
the distance A£ = 2 m of the two points, the bulk con- 
nected correlation function has been computed [10 , and 
it holds 

•41 ^£[<e<e*+A/> (n) - <^> (n) (@W n) ] 

£=1 

= Tr[(9 6) ^ og 2 A V)] = ^ l0S2 Ai h d l °Z2 A*)(13) 

where the summation is performed over the eigenvalues 
f of the CPT map p = \ (V L ®V L +V R ® V R ), and 
where h^(-) is a polynomial function of its argument (in 
this expression a stands for a traceless operator while fi- 
nally N = 2 n+2 is the size of the associated homogeneous 
MERA with periodic boundary conditions). The result 
then follows by noticing that by construction = T> ®T> 
so that the £ can be expressed as products kk' of the 
eigenvalues of V. In particular if as in Eq. ( [l2| ), 6 is 
an eigenvector the Heisenberg adjoint of V we have that 
G © is an eigenvector of the adjoint of p at the eigen- 
value £ = k 2 and thus, 



= cV log 2 Ai = c' (A£) 2l °^ K 



(14) 



which proves the claim (here c' = Tr[(0 0)<r]). 



IV. BOUNDARY CONTRIBUTION TO THE 
GROUND STATE ENERGY 

In the presence of a boundary, the average of extensive 
observables (the ground state energy for example) does 
contain a bulk and a boundary contribution (negligible 
in the thermodynamic limit). In this section, we eval- 
uate the density of the ground-state energy for a local 
Hamiltonian (with interactions among sites at maximum 
distance v) of the form 



n = 



L 



1 



(15) 



where v is the number of sites involved in the model in- 
teraction H. While this problem is easily solved in a 
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level-recursive manner for a MERA structure (in which 
periodic boundary conditions hold) , when explicit condi- 
tions over a defined boundary are involved things become 
slightly more complicated. 

Suppose, for simplicity, that the interaction is again a 
v — 3-body operator, therefore 



(H) 



where 



_(n) 
P3 = 



2 n+2 



Tr[tf 3 -4 n) ] 



E(n) 
Pj,j+lJ+2 



(16) 



We need to build a recursive function which relates this 
average density matrix to the one belonging to the pre- 
vious tensor level p^"^- Of course, the boundaries will 
play some role too in this relation 



P3 



(n) 



1 



2 n+2 



2 n + 1 - 1 



(here V is the average of T>l and Vr). This equation 
shows the contributions of both bulk and edge terms; 
nevertheless, when we approach the thermodynamical 
limit n — > oo, the contribution of the first two terms 
vanishes in every norm, because any density matrix has 
trace norm bounded by one and CPT maps are contrac- 
tive. 

This means that the extensive influence of the bound- 
ary upon the lattice grows slower than the size of the 
system, a physical sounding and known property. To 
quantitatively describe such behavior, we compute the 
(total) energy associated with the block of the first 2 r — 1 
spins near a boundary, say the left one. In our notation 
this corresponds to 



E 



(n) 
1...2 T 



Tr 



2 T -1 

£„(») 

j = l 



(18) 



Exploiting the usual formalism of level-growing CPT 
maps, we can rewrite the sum in (18) as 



2^ Pj,3+1 



i+2 



r-1 

E 

p=0 



(19) 

Now, we can successfully approach the thermodynamical 
limit while keeping r fixed. Recalling that p^ A 1 2 is the 
fixed point of 6l, we obtain 



E 



(th) 

1...2 T 



Tr 



p=0 



(20) 



As expected, the result diverges for r — > oo since the se- 
ries is made of terms growing in trace norm. To explicitly 



estimate how this quantity deviates from its correspond- 
ing value in the bulk as r grows we evaluate the following 
quantity 



AeS t _! = -(2 r - 1) Tr 



Hs-pI 



Tr 



Tr 



tf 3 -^ 2? V?o1C L (p£ >1>2 ) 
p=o 



p=0 



(21) 



Notice that in this case the map V applies to a traceless 
operator, therefore if we decompose the argument in a 
basis of generalized eigenvectors for P, it must have null 
component over the unique state of eigenvalue one. As 
a result the boundary contribution to the ground state 
energy has the form 



AE- 



(th) 



E E 



(22) 



where g(-) is a polynomial in its main argument. Looking 
at Eq. p2| , we notice that the inner sum diverges for any 
eigenvalue kt> of V greater or equal to 1/2, unless the g 
are identically zero for such values of kx>- In general this 
will happen when — JCl(Pa i 2) nas nu ^ component 
over any generalized eigenspace whose k,t> > 1/2. 

Interestingly enough the capability of such deviation to 
diverge is another manifestation that the MERA states of 
Fig. [l] are critical: indeed only the integral of a power-law 
function can diverge, while for an exponential decaying 
correlation function the integral of the corresponding ac- 
tion is always finite. 



V. CONCLUSIONS 

In this paper we exploited the properties of MERA 
to describe boundary critical phenomena. We consid- 
ered the case of a one-dimensional critical system with 
a boundary. To this end we modified the local struc- 
ture of the MERA at the boundary to account for more 
flexibility in its description. 

Besides showing how to compute local observables, we 
achieved two main results. First of all we showed, as pre- 
dicted by boundary conformal field theory, that the crit- 
ical exponents associated to the decay of the one-point 
function (as a function of the distance from the bound- 
ary) is always half of the one of the bulk two-point corre- 
lation function corresponding to the same scaling oper- 
ator. Secondly we compute the boundary corrections to 
the ground state energy and determined its scaling be- 
havior. As in the bulk case, also in the presence of the 
boundary, most of the critical properties are determined 
solely by the eigenvalues of the MERA transfer matrix. 
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A remarkable feature of treating boundary critical phe- 
nomena with MERAs is that it is enough to consider uni- 
form tensor network. This is the result of the scale invari- 
ance of the underlying tensor network which holds also in 
the presence of boundaries. In addition to the practical 
advantage in the numerical simulations, this observation 
further clarifies the properties of the MERA. It is worth 
noticing that such property (that is at the basis of the 
effectiveness of a bounded MERA) is physically equiv- 
alent to the fact that, in boundary critical phenomena, 
the operator content of the bulk is not influenced by the 
boundary [14 , suggesting that maybe the connection be- 
tween MERA and general renormalization group theory 
is deeper than what nowadays understood. 

One dimensional systems display the peculiar feature 
that the boundary can be critical only when also the bulk 



is. This is not the case in higher dimensions [14], where 
we can have a critical boundary in a gapped system, re- 
sulting in a richer scenario for the boundary-bulk phase 
diagram. This richness will reflect in the possibility to 
have different compositions of tensor structure. In this 
paper we considered a matrix product state (at the bor- 
der) connected to a MERA. It is easy to imagine that, to 
describe critical surfaces in a non-critical bulk, different 
compositions of tensor networks are required. 

During the writing of this work a paper by G. Even- 
bly et al. appeared on the archive discussing boundary 
critical phenomena using MERA, see Ref. [24] . 

We acknowledge fruitful discussions with M. Cam- 
postrini, S. Montangero and M. Rizzi, and financial sup- 
port from FIRB-RBID08B3FM and the National Re- 
search Foundation and Ministry of Education Singapore. 
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